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An ultra-large Liquid Argon Time Projection Chamber-based neutrino detector will have the 
uncommon ability to detect atmospheric Vt/Vt events. This paper discusses the most promising 
modes for identifying charged current v-r/V-r, and shows that, with simple kinematic cuts, ~30 v-t-'tT't 
interactions can be isolated in a 100 kt-yr exposure, with greater than 4cr significance. This sample 
is sufficient to perform flux-averaged total cross-section and cross-section shape parameterization 
measurements - the first steps toward using Vt/Vt to search for physics beyond the Standard Model. 



I. INTRODUCTION 



Of all observed Standard Model particles, we have the least direct experimental knowledge of the tau neutrino, 
v^jV-T* ■ A high statistics, dedicated experiment sensitive to VtIVt Charged Current (CC) interactions has never 
' been performed and would add significantly to our understanding of electroweak interactions and be sensitive to 
C^l . certain hard-to-get manifestations of new physics [ij. Indeed, any opportunity to collect a significant sample of 
CC Vt jv-y interactions and simply measure the cross-section and perform basic tests of Standard Model predictions 
' (— I 'I would qualitatively improve our knowledge of the third neutrino weak eigenstate. 

O i' The path to the observation of the Vr was long, not unlike that of the electron neutrino, v^, and the muon neutrino, 
' ] v^. The discovery of the tau-lcpton in 1975 Q led to the assumption of a third neutrino, v^^ the wcak-isospin partner 
of the third charged lepton, r. Since then, indirect information on the has been collected from a wide range of tau 
' decay analyses [3|-|9|. The fact that the is a state orthogonal to and v^.^ for example, was first indirectly revealed 
by the LEP experiments via precision measurements of the Z-boson width [l3|- The direct observation of CC Vt/Vt 
interactions is a very recent (21*^* century) development. The first events were presented by the DONuT experiment 
in 2000 and published in 2001 [ll|- To date, ten CC Vr interactions have been observed, nine by DONuT and one by 
OPERA [H, [l^^ DONuT was a short-baseline, emulsion-based experiment. The Vr/^r^ were produced by a fixed 
QQ ' target 800 GeV proton beam configuration, primarily through the decay of Ds mesons, with the relevant branching 
^\ • ratio B{D^ — !• T'Vr) = (6.6 ± 0.6)% [l^. The Vr and Vr contributed about 3% of the total neutrino flux. 

■ An alternative method for producing/detecting Vt/vt relies on neutrino oscillations. Our current understanding 
QQ , of the neutrino oscillation data [13, 13 indicates that neutrinos produced as v^/v^^s will oscillate mostly into Vr/vrS 
' with an oscillation frequency related to the largest (in magnitude) of the two independent mass squarcd-diffcrences, 
, lAm^gl 2 X 10^'^ eV^. This is true as long as the oscillation length associated with the smallest mass-squared 
difference, Am^j ~ 8 x 10~^ eV^, is much longer than the characteristic baseline of the experiment. The relevant 
mixing angle is consistent with maximal (sin^ 26'23 ^ 1) [13, 113 so a detector placed at one of the oscillation maxima 
and exposed to an originally- beam provides an ideal setup for collecting a large sample of CC Vy jv-r interactions. 
In practice, the design of such an experiment has been demonstrated to be challenging. A large oscillation phase 
demands baselines above (and, preferably, well above) 1000 km, as the tau production threshold requires Vr/VrS. with 
laboratory energies above 3.5 GeV. To date, it has not been possible to produce a beam with an experimentally 
significant neutrino flux at distances beyond about 1000 km [l^[2^. The OPERA experiment [2l|, [l^l, a 1.25 kt 
emulsion-based detector, is aimed at directly observing CC Vt events in a long baseline beam. Unfortunately, the 
L/E ~ 700 km/20 GeV factor does not allow for the collection of a large data sample. As of this writing, one Vr 
candidate event has been observed [l^], with 10.4 events expected after five years of running at design luminosity j2ll |. 

Emulsion detectors provide the strongest resolving power for CC VtIvt interactions. Such detectors isolate the 
events through the observation of a "kink" from the short-lived tau decay. However, automatic scanning of the 
emulsion is a time-intensive process which cannot proceed in real time. An alternative method, with real-time, fast 
event reconstruction using drift chambers was employed by the NOMAD experiment in an oscillation search (23j . 
No Vt interactions were observed as the accessible tvrr? range was outside of what we now know is allowed. Liquid 



* The distinction between Vr and v-r is made now in order to avoid ambiguity later in the paper. 

t This is dwarfed by, for example, the world's growing sample of reconstructed top quark events, which currently consists of well over 
1000 events [H, [ill . 
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Argon Time Projection Chamber (LArTPC) detectors have been proposed for Vt/vt event searches, but have yet to 
be specifically employed for this purpose [ij . Very recently, the ICARUS T-600 neutrino detector started operating in 
Gran Sasso. The ICARUS T-600 is a 600 ton LArTPC exposed, like OPERA, to the CNGS beam [H,!!!]. Although 
the experiment's main goal is to observe disappearance, ICARUS is also equipped to observe a number of CC Vr 
interactions [26j . Unfortunately, the expected event sample is too small to allow one to perform a cross-section 
measurement. In both the NOMAD and ICARUS experiments, the beam direction is used to search for evidence of 
missing transverse momentum consistent with tau production and decay. 

The proposed 20 kt LArTPC at the Deep Underground Science and Engineering Laboratory (DUSEL) opens a 
new opportunity to collect a significant sample of Vr/vr events. Although the detector will be significantly larger 
than OPERA, the combination of longer baseline and significantly lower beam energy leads to anywhere from 0.2-4.0 
expected CC Vr beam-oscillated events/MW-lO^s/kt, depending on the (as-yet-undecided) beam tune-^ [27| . However, 
as a deep underground detector, this experiment will be sensitive to atmospheric neutrino interactions. A natural 
sample of high energy, earth-diameter-as-baseline oscillations, has the drawback of having no clear beam direction. In 
this paper, however, we show that this problem can be overcome so that a significant number of atmospheric v^- jvr 
events can be identified. We present a discussion of the capability of Vt-Wt cross-section measurements an an example 
of the physics potential of this signal. 

The idea of a search for atmospheric Vt jv-r appearance has been pursued before. The 22.5 kt (fiducial) Super- 
Kamiokande detector results disfavor the no Vr/Vr appearance hypothesis at the 2Aa level. However, this 
analysis was hampered by the lack of precision track reconstruction and particle identification inherent to Cherenkov- 
based detectors. A multi-kiloton LArTPC, on the other-hand, will provide well-reconstructed events which can be 
used to isolate the VtI^t interaction signal in a very convincing way, as we present below. We note that the authors 
of [1^ and [sOjUH have also considered the possibility of studying atmospheric Vt/^t using the Ice Cube Deep Core 
Array and a LArTPC (with magnetic calorimeter), respectively. The physics issues involved in the calculation of 
the VtIVt cross-section at hig h energies relevant to Ice Cube Deep Core Array- like experiments have recently been 
discussed in the literature [S^l . 

This paper is organized as follows. In Section [Hi we discuss the "production" of Vt/vt via atmospheric neutrino 
oscillations and review the cross-section and kinematics of CC Vrjvr interactions. In Section IlIIl we present a brief 
overview of LArTPC technology and details of our CC Vt/vt event simulation and the relevant CC and Neutral Current 
(NC) atmospheric neutrino induced backgrounds. In Sections IIVI fVl and IVI[ we discuss the search for CC Vt/Vt 
events via statistical inference in various tau decay modes, and offer some concluding remarks in Section IVlII 

II. ATMOSPHERIC TAU NEUTRINOS: PRODUCTION AND DETECTION 

In order to understand an ultra-large LArTPC's ability to detect it is important to know the properties of 

the atmospheric VtIvt "beam" and the Standard Model expectation for the CC i^T-iV — > tX cross-section, where N 
is a nucleon and X is any hadronic final state. Both issues are discussed in this section. 

A. The Atmospheric Vt/vt Flux 

Most atmospheric neutrinos are a result of pion decays (with a sub- leading kaon component), with the pions 
produced when cosmic rays interact with the Earth's atmosphere several kilometers above sea level. One naively 
expects a I'fi/Vfi to Ve/ve production ratio close to two with virtually zero Vt/vtS^. 

Figure [T] depicts the flux of i/^, F^, i/g, and as a function of the neutrino energy, for cos^zonith = ^1 (left) and 
cos^zonith = (right), where cos6'zcnith is the cosine of the zenith angle^, according to [s^. Figure [2] depicts the 
fluxes as a function of zenith angle for £' = 5 GeV (left) and = 30 GeV (right). The uncertainty on the absolute 
normalization of the atmospheric neutrino fluxes is estimated to be around 20%, while the energy and zenith angle 
dependencies (shapes) are known to within 5% [H, H^l . 



We do not speculate on the possibility of LBNE-beam- based Vt /Vt detection as it is extremely dependent on the largely undecided 

beam parameters (e.g. beam energy spectrum). 
§ A very small Ur/Vr component to the parent atmospheric neutrino flux is expected from, for example, charm production. This addition 

is both very small and of higher energy than the atmospheric neutrinos considered here, and will be neglected henceforth [33Ll34l|. 
1 cosS^enith = corresponds to neutrinos coming from above and cos^^enith = — 1 to neutrinos coming from the antipodal point on the 

Earth. 
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FIG. 1: Atmospheric neutrino flux as a function of the neutrino energy, for cos ^zenith ~ —1 (left) and (right). 
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FIG. 2: Atmospheric neutrino flux as a function of cos^zenith for E = 5 GeV (left) and 30 GeV (right). 



Between production and detection (while the neutrinos traverse a distance L), neutrinos oscillate. The relevant 
vacuum oscillation phases are, 

_ \Aml^\L _ ( lAmf^-l \ A GeV\ / L \ 

AE " 1^2.4 X 10-3 eV2j^ E Al2756 kniy'' ^ ' 

where 12756 km is the maximum diameter of the Earth. $13 can be of order one for atmospheric neutrinos as long 
as L is not much smaller than the Earth's radius. For cos^zenith = 0, i ^ ^JiR^h ^ 440 km (i?® is the Earth's 
radius and h ^ 15 km is the average height above sea- level where neutrinos are produced), and $13 ^ 0.3 for 5 GeV 
neutrinos. $12 is always small unless the neutrino energies are below 1 GeV; $12 < 0.35 for neutrino energies above 
3.5 GeV (the CC Vr/vr production threshold) using Am^2 = '''•6 x 10~^ eV^. In summary, the oscillated Vrlvr 
flux above tau threshold comes from negative cos ^zenith and from the dominant "atmospheric" oscillation frequency, 
proportional to jAm^gj. 

The relevant oscillation probabilities are, ignoring matter effects, 

P^r cos^ 6*13 sin^ 26*23 sin^ $13 --sin2$i3, (2) 
Per ^ SIT? 6*23 sin^ 26li3 sin^ $13 < 0.09 sin^ $13, (3) 

making use of the upper bound on ^13 from Figure [3] depicts Per (left) and P^t (right) as a function of the 

neutrino energy with L = 8000 km (cos^zenith = —0.63) for different values of ^13 and assuming that the neutrino 
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FIG. 3: Pii^e Ur) (left) and P{u^ 
normal neutrino mass hierarchy. 
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(right) with L = 8000 km for sin^ 613 = 0.005 (solid) and 0.01 (dashed), and a 



> 

o 

o 

0) 



10 



1 r 



0.1 



0.01 



0.001 




10 15 20 
Energy (GeV) 



1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 

C0s(6'zenith) 



FIG. 4: Atmospheric Ut and Vt fluxes from oscillations as a function of neutrino energy (left) and cosine of the zenith angle 
cos^zonith (right). In the left panel, the zenith angle is fixed at cos^zenith = —0.4, while in the right panel the energy is fixed 
at i5 = 5 GeV. We assume a normal neutrino mass hierarchy and 613 = 0. 



mass- hierarchy is normal (Am^g > 0). Note that the simplifying approximations, allowing us to write Eqs. ^ and ([3]) 
above, were not made in generating Figure [31 and that matter effects were properly taken into account. The PREA'I 
density profile was used to model the density of the Earth [11]. We safely conclude, further remembering that the 
^ti/i^ii flux is larger than the Ve/ve flux, that the majority (>90%) of the Vrjvr atmospheric neutrino flux at energies 
above a few GeV comes from — > Vr (F^ Vr) oscillations. For large values of 6*13 and a normal (inverted) mass 
hierarchy, there is a small but potentially significant fraction of Vr (F^) from v,, (pe) oscillations. 

Figure m depicts the atmospheric Vr/vr flux at the detector site as a function of energy (left) and zenith angle 
(right), for a fixed zenith angle and energy, respectively. We use the current best fit values of the neutrino oscillation 
parameters - Amf2 = 7.6 x 10"^ eV^, Am^g = 2.4 x lO'^ eV^, sin^ 6*12 = 0.31, and sin^ 6*23 = 0.5 - and 

assume sin^ ^13 = and a normal neutrino mass hierarchy. Prominent oscillatory features are clear in both panels of 
Figure m In summary, Vt/vt-s arrive at the detector mostly from "below" and are overwhelmingly low-energy. Also, 
we expect most of the CC Vt/vt initiated events to occur close to tau production threshold. 



B. The u-rN tX Cross-Section 



In the energy region of interest, CC v-rjvT interactions are well described as if the v^^ were scattering mostly off of 
nucleons. The kinematics of v-rn — > r^p (or Ft-p — r+n) dictate that the minimum neutrino energy in the rest-frame 
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Model 


events 


Vt events 


Total 


NuTeV [44| 


56.9 


24.9 


81.7 


Kretzer et al. [451 


37.7 


17.5 


55.2 


Hagiwara et al. [461 


33.7 


18.5 


52.3 


Paschos et al. [43] 


65.2 


29.6 


94.8 


Nuance [471 


54.1 


23.1 


77.2 



TABLE I: The expected number of CC Vt/vt events/100 kt-yr on an isoscalar target as predicted by various cross-section 
models. The Nuance prediction for interactions on an argon target is also shown. 



of the target nuclcon is 



2m 



N 



= 3.46 GeV, (4) 



where m-r is the tau-leptoii mass and is the nucleon mass, and assuming that the neutron and proton masses are 
the same. 

Similar to CC Vel^e and v^/v^ scattering, CC Vrjvr scattering receives contributions from a variety of processes. 
At low energies, quasi-elastic (QE) scattering Vt-u — > r^p {prV ^ T+n) dominates, while at high energies the deep- 
inelastic scattering (DIS) contribution is largest. At intermediate energies, it is expected that resonance-production 
and other non-perturbative QCD phenomena dominate. The heavy tau mass (m,- = 1.777 GeV) shifts the regions 
where the different contributions dominate towards higher energies. For example, the production of a A-resonance 
requires E^^ > 4 GeV for CC scattering, as opposed to E^^^ > 0.44 GeV in the case of CC v^^ scattering. For 
CC Vr scattering in the region of interest, i?min < ^ 20 GeV, the three different contributions are similar. QE is 
expected to dominate very close to threshold while DIS is dominant above 10 GeV or so [43| . 

We will not discuss the challenges associated with computing the CC Vt/vt cross-section. Instead, we will compare 
different results in the literature in order to illustrate the uncertainty of the situation. While a significant improvement 
is expected in the near future from upcoming experimental data on I'^/l?^ and i^e/Fg-scattering, a large uncertainty 
will remain in the jv^ sector until more data becomes available. 

Table H] contains the expected number of atmospheric CC Vt/vt events/100 kt-yr from various cross-section com- 
putations found in the literature [43l - |46j and the Nuance [l^ prediction. The NuTeV collaboration has extracted Fe 
structure functions from the u-Ye and i7-Fe deep inelastic differential cross-sections using its high energy sign selected 
beam [31 ■ We compare the event rate from this data with those obtained from other cross-section models. The 
calculations done by Paschos et al. [i^ use leading order (LO), whereas Kretzer et al. [i^ and Hagiwara et al. [1^ 
use next to leading order (NLO) structure functions to calculate the total cross-section. The Kretzer et al. model 
takes into account charm production, tau threshold, and target mass effects. The Hagiwara et al. model accounts for 
the effect of polarization on the DIS cross-section, unlike other models. All of these effects are only important in the 
low energy region, however. The significant discrepancy in the high energy region points to the lack of data presently 
available. Figure [5] depicts different predictions for the CC Vr cross-section at neutrino energies below 30 GeV, where 
we expect the vast majority of events to occur. We bring attention to the fact that the different computations in 
the literature find different effective thresholds and that the shape of the rise in the cross-section for energies close 
to threshold also changes from one estimate to the other. We should note that the Monte Carlo model used by 
Super-Kamiokande [2^ has a cross-section very close to that of Nuance, the neutrino event generator used in this 
paper and discussed later. 

In an attempt to quantify these differences, we simply parameterize the cross-section using 



1 - exp [-aVE^]) X IQ-^sfffL, (5) 



E V V ' ^jj -^- 

where the parameter h fixes the production threshold value in the lab frame and the parameter a fixes the rate of 
increase near the threshold. The energy E and parameter h are in GeV. We find that this parametrization provides 
an excellent fit to the curves shown in Figure [SJ The best fit values for these parameters are tabulated in Table [Hi 
Later, we calculate the experimental sensitivity to these parameters, assuming a jvr measurement consistent with 
our Monte Carlo expectation. 
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FIG. 5: The CC v-r cross section for an isoscalar target calculated according to different models: Paschos et al. (long-dashod), 
Kretzer et al. (solid) and Hagiwara et al. (short-dashed). The Nuance (dotted) cross section prediction for an argon target is 
also shown. 



Model 


a 


h 


NuTeV [44] 


0.085 


3.9 


Kretzer et al. [451 


0.080 


5.6 


Hagiwara et al. [46] 


0.089 


6.7 


Paschos et al. [43] 


0.105 


4.1 


Nuance [47] 


0.089 


4.7 



TABLE 11: The best fit parametrization for the various CC Vr cross-section models. See Equation [5] for definition of o and h 



III. EXPERIMENTAL SETUP AND EVENT SIMULATION 
A. Liquid Argon Time Projection Chamber Technology 

The LAR20 detector provides a typical design for an ultra-large LArTPC [H]. We assume a 100 kt-yr exposure for 
our event rate calculations, a modest estimate with the proposed 17 kt fiducial volume design and >10 years of data 
taking. Three wire planes, each with 3 mm wire spacing and 2 MHz sampling, will be instrumented. This granularity 
is insufficient to directly reconstruct the tau decay kink and so the analysis presented here relies on indirect methods 
of isolating CC Vt/vt interactions. The TPC's active volume will be 15 x 15 x 34 m'^, featuring multiple field cages 
with voltages of 500 V/cm (corresponding to a drift velocity of 1.5 mm//is). Modest shielding is required for such 
a large detector so that the atmospheric neutrino events are not masked by and/or confused with cosmic ray tracks 
and interactions. LAR20 is proposed for the 800 foot level of DUSEL, which provides more than enough shielding for 
this analysis. Although a magnetized LArTPC has been demonstrated to work ji^, the current LAR20 design does 
not include this feature. It should be noted that a magnetic field would greatly enhance the analysis described here, 
allowing for possible and Vr separation and finer energy resolution, among other things. 

The spatial resolution for reconstructed tracks in the detector is at the millimeter scale in the drift and wire 
directions. A particle's energy can be reconstructed by adding up all of the charge collected along a stopping track. 
In the case that the particle is identified with some confidence, measuring the distance of range-out and/or multiple- 
scattering (sol can also be used to reconstruct the energy of the particle. These techniques can be used simultaneously 
to measure the energy more precisely. Although the actual energy resolution is dependent on the particle's energy and 
identity, the energy resolution for LArTPCs is usually quoted at the few percent level [5lj. An example of data-based 
LArTPC neutrino event reconstruction can be found in Ref. [s^ . In this analysis, we consider energy resolutions of 
0% and 15%, reconstructed zenith angle resolutions of 0° and 10°, and 100% charged particle identification efficiency. 

Although precise energy/angular resolution is not particularly vital to this measurement as demonstrated below, 
the identification of charged and neutral pious is. Especially at >1 GeV energies, separating a charged or neutral 
pion from a proton/electron/muon/kaon, is nearly 100% efficient with LArTPCs. Particle identification proceeds in 
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Liquid argon 


Nucleon binding energy = 29.5 MeV 
l^ermi momentum (pj — 242 MeV 
Fermi momentum (n) = 259 MeV 
Density = 1.396 g/cm^ 


Nuclear density 
Pir) = Twfe 


c = 3.53 fm 
z = 0.542 fm 



TABLE IIL Some relevant parameters in the simulation of neutrino interactions on (liquid) argon. 



several ways. First, there is the combination of energy deposition per unit distance (dE/dx) and range. dE/dx itself 
can distinguish highly ionizing particles such as protons and kaons from minimum ionizing (^2.1 MeV/cm in liquid 
argon) particles such as charged pious and muons with close to 100% efficiency [53j . High energy electrons and gammas 
(usually from neutral pion decay) are identified confidently as they lose most of their energy via Bremsstrahlung and 
e"^/e^ pair production, creating a well defined electromagnetic shower. The experimentalist can also take advantage 
of the gamma's 18 cm conversion length in liquid argon, often resulting in a discernible gap between the interaction 
vertex (gamma creation point) and the beginning of the shower. 

Pions are separated from muons via hadronic multiple scattering (with hadronic interaction length of 84 cm), 
nuclear capture, and decay products. A particle that travels longer than a few hadronic interaction lengths without 
a secondary interaction can be identified as a muon. In the case that a negatively charged particle does not decay 
in flight, a tt" (^~) will capture on argon 100% (76%) of the time. Negatively charged pions and muons can be 
differentiated in this way as the nuclear capture products are significantly different in each case. The decay chains 
(tt — > /i — > e) and (/i — > e) can also be used to separate pions from muons. 

The majority of atmospheric Vr/vT interactions are high Q^, deep-inelastic scattering events usually with large 
multiplicity and vertex activity. Such events can be difficult to fully reconstruct for Cherenkov-based experiments and 
detectors with weak spatial resolution. Disentangling long, energetic tracks (possibly overlapping rings in the case 
of a Cherenkov-based detector) is relatively simple for LArTPCs featuring three-dimensional imaging and millimeter 
resolution in a homogeneous and fully active volume. Even in the case of multiple 7r° production and subsequent decays 
(7r° — >■ 77), gamma pair matching is straightforward with the three-dimensional imaging capabilities of LArTPCs and 
the help of a reconstructed invariant tt" mass. 



B. Neutrino Event Generation 

The Nuance (version 3) neutrino event generator [47| has been employed to simulate all-flavor atmospheric neutrino 
interactions on argon. The Nuance source code, originally created for simulating the interactions of atmospheric 
neutrinos with water, has been modified slightly in order to properly simulate the neutrino- argon interaction and the 
propagation of the resulting hadrons through the argon nucleus. The reader is referred to the Nuance documentation 
for detailed information on the various cross-section and nuclear process models used in the program. Modeling 
the energy and angular distributions of the products of intra-nuclear collisions is difficult as there are many types 
of interactions and little data. In Nuance, hadrons are stepped through an argon nucleus with measured radially- 
dependent density distribution and Fermi momentum. Some relevant argon-specific parameters are listed in Tablelllll 
The pion- nucleon and nucleon- nucleon cross-sections and angular distributions are largely based on HERA data |54| . 
The understanding of the cross-sections and angular distributions of relevant processes such as pion absorption, charge 
exchange, and elastic and inelastic scattering within the argon nucleus will be greatly improved by the ArgoNeuT [55| . 



MicroBooNE [56|, and ICARUS T-600 l2J] experiments in the near future. These experiments will also improve our 
knowledge of the neutrino-on-argon cross-sections themselves. Similarly, MINERi^A [s^] will greatly enhance our 
knowledge of cross-sections and intra-nuclear interactions using multiple nuclear targets in the near future. The 



authors also look forward to a full intra-nuclear cascade simulation from the GENIE collaboration [58[, featuring the 
simulation of "all (intra-nuclear) reactions on all nuclei" . 

Of special relevance to this paper, Nuance decays tau-leptons with the TAUOLA (version 2.6) package [11] • In CC 
quasi-elastic and deep- inelastic Vt I^t interactions, the polarization of the tau is calculated using the appropriate form 
factors [60j . The tau is considered com plet ely polarized in resonant interactions, where the form factors are usually 
ambiguous and/or difficult to calculate |47| . 
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IV. HADRONIC TAU DECAYS, INCLUSIVE 

Identifying a Vt /T^t on an event-by-event basis will be extremely difficult for a detector unable to resolve the kink 
from the tau decay, occurring only a fraction of a millimeter from the event vertex for the most relevant atmospheric 
neutrino energies. However, statistically inferring the presence of CC Vt/vt interactions in a capable detector is 
possible by analyzing event kinematics. 

The tau decays to one or more hadrons and a neutrino about 65% of the time, with the branching ratios of the 
various channels well known. For convenience, we tabulate the relevant decay modes in Table IIVI Looking for 
exclusive hadronic final states may assist in confirming the Vt/i^t appearance hypothesis. However, an exclusive 
analysis (requiring a specific pion final state) is complicated by an increased dependence on the simulation of the 
relatively uncertain final state interaction processes. We briefly consider the exclusive channels in the following 
section. In this section, we consider both inclusive pionic decays of u-r and Vt together. Note that there is no reason 
that kaonic decays, compromising only about 3% of the total hadronic channel, cannot be included in this sample. 
However, we do not consider them for simplicity and because of their negligible effect on the analysis. 

The largest background to Vt/vt detection via hadronic-channel-inclusive pion decay kinematics is the atmospheric 
neutrino NC interaction involving the excitation(s) and then subsequent decay(s) of a nuclear resonance to a pion 
and a nucleon (with a small pion contribution from elsewhere, such as NC coherent production). CC v^^jv^ and 
Ve/ve events are considered a negligible background for Vr /^r detection via the hadronic channels as the efficiency 
of GeV-scale electron and muon identification is near 100% for LArTPC-based neutrino detectors [s^l. Separation of 
the CC and NC channels is critical to the detector's search for Ve/ve appearance in the accelerator-based part of the 
experiment. 

Pions originating from neutrino-induced nuclear resonance decays (and elsewhere) and pions from tau decay form 
considerably different kinematic distributions. Variables such as visible energy, average energy per pion, angle between 
the highest energy pions, angle between the highest energy pion and the (reconstructed) initial neutrino angle, energy 
of the highest energy pion, invariant mass of the pionic system, and more are all useful, correlated, observables in 
differentiating pions originating from Vt/vt events and pions from the NC background. 

For a 100 kt-yr exposure, 77.2 CC i^r+i^r events are expected, based on the previously discussed oscillated atmo- 
spheric neutrino flux and the Nuance cross section predictions. Of these, 46.3 events feature at least one charged pion 
from tau decay. CC jv-T events involving a charged Icpton from the tau decay are not considered here, even if a 
non-tau-dccay-pion escapes the nucleus. Visible energy (i?vis), cos(reconstructed zenith angle) (cos(6'zcnith(ioc.))) ^'^d 
energy of the highest energy neutral or charged pion (i?^7r) have been found to be the most useful variables in differ- 
entiating signal from background. These variables also happen to be among the simplest to measure experimentally. 
We initially consider events with at least one visible charged pion. Note that the pions in the CC Vt^Vt sample are 
not all necessarily from the tau decay (e.g. they might be resonant-induced) although a tau decay to at least one 
charged pion is required to enter the sample. 

Events with cos(0zonith(rec.))<~O.2 and one or more charged pions enter the i^^+Ft- candidate sample. Reconstructed 
zenith angle is defined assuming full knowledge of all final-state and tau-daughter non- neutrino/neutron particles. 
The reconstructed zenith angle cut is applied to reduce the downward-going NC background as few Vrl^r events are 
expected in this region. The cut is stricter than the usual cos(0zenith(rec.))<O-O cut in order to compensate for the 
limitations of reconstructed angle resolution in differentiating upward- from downward-going events. 

With a 3.5 GeV threshold, CC Vr/vr interactions generally occur at higher energies than the other possible atmo- 
spheric neutrino interactions. Visible energy is defined as the sum of the energy of all primary particles and tau decay 
products minus the energy contribution of any and all neutrinos and neutrons. Figure [5] shows the visible energy of our 
sample, separated into NC, CC Vr+Vr, and total after simulating atmospheric neutrinos with energy 3<-B,y<30 GcV. 
Requiring i?vis>6.0 GeV is found to be effective in isolating the Vr/vr events from the NC background. The visible 
energy cut was chosen to improve signal-to-background, although a thorough optimization of this cut has not been 
performed. Ref. [30| also found a 6-7 GeV visible energy cut useful for isolating atmospheric v-t jV-Y events. 

Figure [7] shows E^-^ for each sample after the visible energy cut. After further requiring i?^^>1.5 GeV, we are 
left with 73.4 total events over a background of 44.8 events, a 4.3cr excess (with Poisson statistics only) assuming 
experimental data that is consistent with the Monte Carlo expectation. This excess has been found to be largely 
insensitive to the effects of even modest energy and reconstructed angle resolutions. Applying a Gaussian smear (with 
a 15% sigma) to each event's i?vis and Ea^t^ and a Gaussian smear (with a 10° sigma) to cos(6'2(,nith(rec.)) (all smears 
together), left the excess at 3.9(t. Similarly, varying each cut individually by ±15%, while leaving the other two cuts 
at their nominal values, dropped the excess down to 4.2(t at worst. 

The background event rates quoted above rely on the understanding of the NC cross sections in the ~few to 30 GeV 
energy range and final state interactions. Fortunately, experiments such as those mentioned in Section IlII BI will soon 
improve our knowledge of these processes. The Nuance event generator represents our best estimate of the NC 
background and wc employ it with the understanding that the "real" background (and signal) may be substantially 
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TABLE IV: Branching ratios of selected r decay modes [16 
mediated by a resonance, '/i^''' is either a pion or a kaon. 



^All but a small fraction (0.3%) of this branching ratio is 



signal (CC i^r +T't) 
background (NC) 
total 




10 12 14 16 18 20 

(GeV) 



FIG. 6: The visible energy for NC interactions with at least one final state charged pion and CC Ut+Vt interactions with the 
tau decaying to at least one charged pion. Events with i5vis>6.0 GeV enter the final sample. 



different from that which is predicted. The apparent flexibility of the resolution/cuts described above demonstrates 
that the sensitivity quoted will only minimally depend on the modeling of NC cross sections and final state interactions. 
Moreover, a strong excess remains even in the case that Nuance severely under-predicts the background. Arbitrarily 
increasing the background that passes all cuts by 50% (while leaving the signal at its nominal level) leaves the excess 
at 3.5(7. 

This analysis would benefit from a neural network and/or likelihood analysis involving the many correlated kinematic 
variables useful for discriminating CC Vt/^t from the NC background. Although a complicated multi-layer approach to 
inferring atmospheric Vt jV-,- appearance could be utilized to improve sensitivity, we have shown that simple observables 
(and a finely-grained and efficient detector with particle identification capabilities) are all that is necessary for this 
measurement. 



V. HADRONIC TAU DECAYS, EXCLUSIVE 



The inclusive hadronic decays are broken up into their constituent exclusive channels in Table |Vl The three cuts 
discussed above (i?vis>6.0 GeV, i?-t-,r>l-5 GeV, and cos(6'5,(,„ith(rec.))<— 0.2) are used to arrive at the expected number 
of events. The categorization of the inclusive channels cannot simply be made according to the expected tau decay 
products and branching ratios as at least one non-tau-decay-pion appears in 68% of the CC Vt jvr events simulated. 
For example, requiring the invariant mass of a charged and neutral pion pair to match the p meson peak, in an effort 
to isolate the most likely decay channel (r n^n^Vr, with a branching ratio of 25.4%), is complicated by these 
non-tau-decay-pions. Figure [8] shows the invariant mass of the highest energy charged and neutral pion in events with 
>0 charged pions, 1 neutral pion, £^vis>6.0 GeV, and cos(0zcnith(rcc.))<~O.2. Although a clear p peak is visible, signal 
events often do not reconstruct well to the p mass. Furthermore, the NC background contains a significant amount 
of p resonances, a consequence of vector meson dominance. Requiring the highest energy charged and neutral pion 
invariant mass to be greater than 500 MeV along with the same three cuts on events with a single neutral pion and 
one or more charged pions only negligibly increases the sensitivity. 

Although the actual sensitivity to Vr/vr appearance cannot be improved via an exclusive analysis, such an analysis 
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Et. (GeV) 

FIG. 7: The energy of the highest energy pion, neutral or charged, for NC interactions with at least one final state charged pion 
and CC Vt+Vt interactions with the tau decaying to at least one charged pion. Non-tau-decay-pions (e.g. resonant- induced) 
are also considered in the CC Vt+Vt sample. Events with _E^^>1.5 GeV enter the final sample. 



Exclusive channel (vr's) 


Total/Background 


Excess significance (cr) 


1 charged, neutral 


2.3/1.6 


0.6 


1 charged, 1 neutral 


5.9/3.2 


1.5 


1 charged, >1 neutral 


6.9/3.5 


1.8 


>1 charged, neutral 


11.5/7.6 


1.4 


>1 charged, 1 neutral 


21.5/14.1 


2.0 


>1 charged, >1 neutral 


25.3/14.8 


2.7 


Total 


73.4/44.8 


4.3 



TABLE V: The expected number of CC u-r-'rVT and NC background events/100 kt-yr after cuts for hadronic tau decay channels 
categorized by number of charged and neutral pions. 

can help to confirm that a possible excess is in fact from Vr/vr interactions rather than non-standard neutrino 
interactions. That is, seeing the expected excess in each charged/neutral pion category consistent with tau decay 
(and the expected non-tau-decay-pion contribution) can rule out non-standard neutrino interactions as the source of 
an excess. 

It is worth noting that atmospheric CC VtIvt events with a kaon in the final state (i.e. r — > Kv) should not be 
considered a background for the proton decay channel p — > K^v. An argon-bound-proton's decay kaon is expected to 
have a momentum of less than 500 MeV/c as it exits the nucleus, taking into account the effects of kaon rescattering 
(with or without an argon spectral function) [6]| . Conservatively requiring the kaon's momentum to be <5Q0 MeV/c, 
less than 0.05 CC Vrjvr ij Kv) cvcnts/100 kt-yr are expected. 

VI. LEPTONIC TAU DECAYS 

The tau decays to a charged lepton and two neutrinos with a branching ratio of about 35%. Accelerator-based 
Vt/vt appearance experiments without the ability to see the decay kink in such events can infer the presence of a 
VtI^t interaction by searching for missing transverse energy (carried away by the two unseen neutrinos). This missing 
transverse energy is absent for background CC ff^/yfi/ve/T^e events which don't feature final state neutrinos. However, 
an atmospheric Vr appearance search is not afforded the luxury of a known beam direction and a missing transverse 
energy search is therefore difficult. 

Compared to the inclusive hadronic modes discussed above, the CC Vr/v-r leptonic channels have a higher back- 
ground and smaller signal. The leptonic channel's background (CC t'^/i^^/z^e/^^e) cross-section is almost three times 
higher than its NC counterpart and the signal tau decays to a charged lepton about half as often as it does to one or 
more hadrons. Furthermore, there are fewer cuts available to separate background from signal. The powerful visible 
energy cut used for the inclusive hadronic final states above is rendered mostly useless for the leptonic case as the 
CC VtIvt^ unseen neutrino daughters push the event's visible energy down, overlapping more with the background 
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FIG. 8: The invariant mass of the highest energy charged and neutral pion for NC interactions and CC Pt+T't interactions 
with one neutral pion and at least one charged pion. Non-tau-decay-pions are also considered in the CC Ut+Vt sample. 
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FIG. 9: The visible energy distribution of charged leptons from signal CC Vt+Vt and background CC v^+V ^^Ve+Vc. events. 



CC y^./V^/ue/Ve visible energy. Requiring cos(6'zonith(rcc.))<~0.2 and i?vis>3.0 GeV, we find 20 CC Vt+Vt events 
with an expected background of 610 CC i^^+T/^+i/e+J^e events (see Figure [S]). Cuts involving the charged lepton's 
transverse momentum and direction, both with respect to the reconstructed neutrino direction, have been found to 
be largely ineffective. We find the leptonic channel an unlikely source for atmospheric Vr/vr investigation. 



VII. CONCLUSIONS 



The Vt/T^t is the least well understood observed Standard Model particle with only ten events ever identified. 
Simply measuring the cross-section and testing the Standard Model prediction, after obtaining a sizable sample of 
events, will improve our understanding of this elusive particle. The prospect of measuring the atmospheric neutrino 
mixing parameters from a j/^ — > i^,- {V ^ — > Ft-) appearance experiment, however imprecisely, rather than & v^^ 
(T^fj. Vfj) disappearance experiment would also be a strong corroboration of the three neutrino mixing model in 
general. 

Along with offering an attractive detector for high- and low-energy accelerator-based neutrino oscillation experi- 
ments, atmospheric v^/v^ disappearance measurements, proton decay searches, and sensitivity to supernova burst 
and diffuse neutrinos, among other things, a kiloton-scale LArTPC will also be capable of observing a significant 
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FIG. 10: The CC i^t cross-section parametrization sensitivity (see Equation [Sjl , assuming a measurement consistent with the 
Nuance prediction for neutrino-on-argon scattering. The lines denote the 1-a (statistical-error-only) allowed region. The other 
points denote models by Hagiwara et al. (H), Kretzer et al. (K), Paschos et al. (P) and the NuTeV collaboration (NT) for 
neutrino-on-isoscalar-target scattering. Note that the confidence level drawn is optimistic as it assumes full knowledge of the 
Vr'-i^T ratio. 




0.25 



FIG. 11: The flux-averaged total CC v-t+Vt cross-section (6<iJvis<30 GeV), assuming a measurement consistent with the 
Nuance neutrino-on-argon prediction (with 1-a statistical error). We show the cross-section for neutrino-on-isoscalar-target 
from other models for reference. Letters adjoining the dots refer to the same models as in Figure 1101 



number of atmospheric Vt/i^t appearance events. Contrary to naive expectations, we find that the most promis- 
ing way of detecting atmospheric v^jv-rS is through the study of the tau hadronic decay modes as it is difficult to 
separate tau-decay leptons and leptons coming from CC Ve/ve and Vy^jv^ interactions. After simple cuts on visible 
energy, reconstructed zenith angle, and energy of the highest energy pion in atmospheric neutrino events with no 
charged lepton, 28.5 tau neutrino events over a background of 44.8 events are expected from a 100 kt-yr exposure, 
corresponding to a 4.3(7 excess signal with Poisson statistics only. Although the detection of such events might be 
considered "indirect" , the expected signal-to-background ratio over most of the kinematic range (after cuts) is >0.5 
and even > 1 at high energies (see Figure [7|) . Furthermore, by studying the exclusive decay channels categorized by 
number of charged/neutrino pions, the possibility of an excess in the inclusive analysis being a result of some form of 
non-standard interaction, rather than u-r /T't production, can be disfavored. 

Given a measurement consistent with the simulated sample, the expected statistical-error-only sensitivities to the 
CC Vt cross-section parametrization and the flux-averaged total CC Vr+Vr cross-section (6<i?vis<30 GeV) from a 
100 kt-yr exposure are shown in Figures [TU] and [TT] along with the predictions of various cross section models. Such 
measurements would provide basic tests of the Standard Model and vastly improve our knowledge of the elusive third 
neutrino weak eigenstate. 
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